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We apply the density-functional technique to determine the lattice constant of GaAs supercells 
containing MnGa, Mni nt , and AsGa impurities, and use a linear interpolation to describe the de- 
pendence of the lattice constant a of Gai-ajMn^As on the concentrations of these impurities. The 
results of the supercell calculations confirm that MnGa does not contribute to the lattice expansion. 
The increase of a is due to both Mnint and AsGa, that are both created in the as-grown (Ga,Mn)As 
in proportion to x, and that are most probably present in a remarkable amount also in the best 
annealed materials. 

PACS numbers: 71.15.Ap, 71.20.Nr, 71.55.Eq, 75.50.Pp 



I. INTRODUCTION 

Diluted magnetic semiconductors (DMS), represented 
by (Ga,Mn)As mixed crystals, are particularly interest- 
ing materials because of their hole-mediated ferromag- 
netism This, combined with the semiconducting 

behavior, makes the DMS attractive for applications in 
spin electronics, and stimulated extensive theoretical and 
experimental studies of these materials in the last years. 

Recently, the compositional dependence of the lattice 
constant of the (Ga,Mn)As mixed crystals has been inves- 
tigated with respect to Mn incorporation in the lattice. It 
is known for long that the lattice constant of (Ga,Mn)As 
increases with increasing content of Mn Q. Originally, 
it was attributed to the substitution of Mn for Ga in the 
cation sublattice and it was assumed that the lattice con- 
stant of (Ga,Mn)As extrapolates to the lattice constant 
of the zinc-blende MnAs crystal that is larger than the 
lattice constant of GaAs 4]. This assumption is, how- 
ever, in a contradiction with the simple estimate based 
on the atomic radii of the constituent atoms because Mn 
atom is smaller (0.117 nm) than Ga atom (0.125 nm) [3. 
Also the density-functional study of ev-MnAs showed jfj 
that the lattice constant of MnAs with a perfect zinc- 
blende structure is smaller then the lattice constant of 
GaAs. 

Most recently, the density-functional method of deter- 
mining the lattice constant by minimizing the total en- 
ergy has been applied also to the (Ga,Mn)As. The mixed 
crystals with realistic concentrations of Mn up to 10 per- 
cent were treated by means of the LMTO method within 
the coherent-potential approximation (CPA) 0. It was 
found that the substitutional Mn (MnQ a ) has only a neg- 
ligible effect on the lattice constant in this concentration 
range. It was deduced that the observed increase of the 
lattice parameter is a secondary effect caused by creation 
of an increasing number of compensating donors, propor- 
tionally to the content of Mn. In particular, the As anti- 
site defects (Asc a ) and Mn atoms in the interstitial posi- 
tions (Mni n t) were shown to cause a remarkable increase 
of the lattice constant of (Ga,Mn)As. Quantitatively, 
the dependence of the lattice parameter a(x s ,Xi,y) on 



the partial concentrations of MnQ a Ga, Mn; nt , and Asc a 
was parametrized in the form 

a(x s ,Xi,y) =a a + 0.002^ + 0.105^ + 0.069y (nm), (1) 

where a a is the lattice constant of the pure GaAs. 

In strongly compensated materials it is expected that 
the number of Mn; nt and Asc a increases proportionally 
to the total concentration x = x s + Xi of Mn 8]. The 
dependence of a on i can be simply estimated from Eq. 
(1) in the limiting case of the complete compensation, 
i.e. for x s = 2xi + 2y. In this case, we obtain a linear 
dependence a(x) k, a a + 0.035a; (nm). The coefficient 
will be smaller for partial compensation which is in a 
reasonable agreement with the experimental value 0.032 
nm Q. 

The theoretical prediction that the increase of the lat- 
tice parameter a with Mn concentration is, at least partly, 
connected with the presence of the interstitial Mn, was 
also confirmed by recent measurements 0, 0, 111 m 
which the lattice constants for as-grown and annealed 
materials were compared. During the the post-growth 
thermal treatment, the interstitial Mn atoms diffuse out 
of the material and the lattice constant decreases accord- 
ingly. These measurements, however, differ in several re- 
spects from the calculations. First of all, the calculated 
increase of the lattice constant due to Mn; nt according to 
Eq. (1) seems overestimated by a factor of 2. Second, 
even the best materials with minimum compensation (i.e. 
with all intcrstitials removed) have their lattice constant 
larger than GaAs and increasing with the Mn content. It 
is interpretted that - after all - the Mn substitution itself 
can expand the lattice. 

That is why we re-examine our original coherent- 
potential study by using the full-potential linearized- 
augmented-plane-wave method (FPLAPW 0) that 
overcomes some simplifications involved in the LMTO- 
CPA study (atomic-sphere approximation, unrelaxed lat- 
tice, etc.). We apply it to the supercells of GaAs with 
MnQa, Mn; nt , and Asca impurities. 
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Figure 1: Calculated lattice constant as a function of the 
concentration of the impurities: (a) Mn atoms in the substi- 
tutional positions (circles), (b) Mn atoms in the interstitial 
positions (triangles), (c) As antisite defects (boxes). The full 
triangle corresponds to a relaxed geometry. 



II. RESULTS 

We use tetragonal supercells formed by 8 and 16 molec- 
ular units of GaAs. These supercells, containing a single 
defect (MnQa, Mn illt , or AsQa), represent materials with 
12.5 and 6.25 atomic percent of the impurities, respec- 
tively. The lattice constant was determined by minimiz- 
ing the total energy E tot of the supercell with respect to 
a, either with or without lattice relaxation. In practice, 
the dependence of E to t on a was approximated by a cubic 
polynomial fitted in approx. 10 points. 

The results are summarized in Fig. 1. First of all, the 
calculations confirm that the substitutional MiiQa atoms 
do not contribute to the expansion of the lattice. The 
dependence of the lattice constant on x s is weak and, in 
contrast to the CPA results, decreasing. Assuming that 
the lattice constant of (Ga,Mn)As follows the Vegard's 
law, i.e. that the function a(x s ) is linear in the entire 
concentration range, we obtain the linear coefficient « - 
0.005 nm. This value is in a reasonable agreement with 
the interpolation between the lattice constants obtained 
for GaAs and a-MnAs 0. 

The interstitial Mn atoms were considered in both 
tetrahedral positions, T(As4) and T(Ga4). Although the 
local relaxation of the lattice is different for these two 
defects [l3j], the effect on the lattice expansion is almost 
the same for both. That is why only the data for T(As4) 
are presented in Fig. 1. Without relaxation (empty tri- 
angles), the lattice expands significantly in presence of 
Mriint and the lattice constant increases almost linearly 
with Xj. If the relaxation around Mn^t is taken into ac- 
count, the expansion is slightly weaker, as indicated by a 
full triangle in Fig. 1. It is important to notice that the 
linear coefficient is « 0.048 nm now, approximately one 
half of the value obtained from the CPA calculations. 

The contribution of the Asc a antisite defects to the lat- 
tice expansion is very similar to Mn; nt . Assuming that 



the contributions of various defects to the lattice expan- 
sion are additive, we summarize the above results into a 
simple linear formula, 

a(x s ,Xi,y) = a a - 0.005x s + 0.048xi + 0.046j/ (nm). (2) 

III. DISCUSSION 

Eq. (2) is an full-potential counterpart of Eq. (1) 
that was obtained from the TB-LMTO-CPA study 0. 
Because of several simplifications involved in the CPA 
calculations (effective medium, mimimum basis, atomic- 
sphere approximation, etc.) Eq. (2) should be considered 
more reliable. At the same time, however, the sensitivity 
of the parameters of Eqs. (1,2) to the method of calcu- 
lations indicates that also the coefficients in Eq. (2) rep- 
resents only a rough quantitative estimate for the com- 
positional dependence of the lattice constant and should 
be used with a caution. 

The most important features of Eq. (2) is that it con- 
firms the recent results |6j, []] that the lattice expansion 
of (Ga,Mn)As can not be attributed to the substitutional 
Mn. It is related to the presence of the interstitial Mn and 
Asca antisite defects as well. The contribution of these 
two intrinsic defects to the expansion is characterized by 
a linear coefficient w 0.05 nm instead of much larger val- 
ues resulting from the CPA calculations. The present 
results are much closer to the experiment |y, EU uM ■ 

To obtain a more detailed picture of the lattice ex- 
pansion of real materials and to compare the role of 
Mnj n t and Asc a we need to estimate the partial con- 
centrations of MnQa, Mnj nt , and Asca- To do this, we 
start with the concentration dependent formation ener- 
gies E(Mn Ga ), £(Mn in t), and E(As Ga ) H Q and as- 
sume that the corresponding partial concentrations in 
the as-grown (Ga,Mn)As can be approximately calcu- 
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Figure 2: Concentration dependent lattice constant of 
(Ga,Mn)As according to Eqs. (2] and (4). (a) as-grown mate- 
rial with MnGa, Mnint, and AsGa (triangles), (b) mixed crys- 
tal with all Mnint removed (boxes), (c) material with the Xi 
reduced to one third (circles). 
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lated from the dynamical-equilibrium conditions |8|, Il5| 

£(Mn Ga ) = £(Mn int ), £(As Ga ) = 0. (3) 

Using for simplicity the linearized form of the formation 
energies we obtain that, except for the lowest values of x, 
the partial concentrations of substitutional and intersti- 
tial Mn, and also the concentration of the Asca antisite 
defects increase proportionally to x, 

x s ps 0.85a;, x t ps 0.15a;, y ps 0.11a;. (4) 

Combining Eq. (2) with the estimate Eq. (4), we arrive 
to the theoretical model of the concentration dependence 
of the lattice constant in the as-grown (Ga,Mn)As. It is 
shown by solid line in Fig. 2. The dotted and dashed 
lines in Fig. 2 correspond to the annealed materials with 
a reduced amount of the interstitial Mn, but with the 
unchanged concentration of the more stable MnQ a and 
Asca- If all Mni„t are removed, the lattice constant is 
still an increasing function of x, but its slope is much 
smaller. An intermediate case with the number of the 
Mn interstitials reduced to one third is considered as a 
realistic example. Most probably, the observed increase 



of the lattice constant in the annealed materials 0, [ToL 
ITll | indicates that a large number of the compensating 
antisite defects (AsQ a as well as some residual Mnj nt ) 
are present even in the best annealed samples with the 
minimum compensation. 

To summarize, the full-potential supercell calcula- 
tions confirmed qualitatively the recent results of the 
CPA studies Q of the concentration dependence of 
(Ga,Mn)As mixed crystals. The present results are 
closer to the experiment. They show the negligible 
influence of the substitutional Mn on the lattice constant 
and the role of the intrinsic compensating donors in the 
lattice expansion. 
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